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Abstract 

We analyze a simple approximation scheme based on the Morita-approach for 
the example of the mean field random field Ising model where it is claimed to be 
exact in some of the physics literature. We show that the approximation scheme is 
flawed, but it provides a set of equations whose metastable solutions surprisingly 
yield the correct solution of the model. We explain how the same equations appear 
in a different way as rigorous consistency equations. We clarify the relation between 
the validity of their solutions and the almost surely discontinuous behavior of the 
single-site conditional probabilities. 
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1 Introduction 

The Morita-approach or equilibrium ensemble approach to systems with quenched dis- 
order goes back to |18) . A fair and clear review from a theoretical physicist's point of 
view containing a quick outline of the theory and various interesting recent applications 
is given by Kiihn in ^1] (see also |13j). 

The central idea in the Morita-approach is: Look at the joint measure governing 
the distribution of the quenched degrees of freedom and the dynamical variables, rather 
than directly trying to describe the quenched measure for the dynamical variables for 
fixed realization of the disorder. Ideally one would then like to write this joint measure 
as a formal equilibrium average over the joint variables in terms of a joint Hamiltonian. 
This joint Hamiltonian would then be the sum of the original one and another "disorder- 
Hamiltonian" depending only on the quenched degrees of freedom. The resulting model 
possesses full spatial symmetries and one might hope that it is amenable to techniques 
known from systems without disorder. 
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Mathematically there are problems with this idea. In finite volume this "disorder- 
Hamiltonian" can in principle be chosen in such a way that the resulting joint distri- 
bution coincides with that of the true model. For lattice systems in the infinite volume 
this is however a serious problem. In fact, for many models an absolutely summable 
joint Hamiltonian does not exist, and the joint measures in the infinite volume are 
non-Gibbsian measures. The appearance of a non-Gibbsian joint distribution was first 
discovered in the example of the Grising model in [3] and studied in a general context in 
|H1 EI] • See also the discussion in EES] • A well- understood example for this are in 
particular the joint measures of the random field Ising model [3] in more than 3 dimen- 
sions at low temperature and small disorder. They provide an illuminating example of 
strong non-Gibbsian pathologies. In fact, their conditional probabilities are shown to 
be discontinuous functions of the conditionings, for a set of conditionings with (joint) 
measure one ^Uj. This means that the measure is not even "almost surely Gibbs" (in 
the sense of |17j). This pathology even causes the usual Gibbs variational principle 
to fail Close analogies to this behavior on the lattice can be already found in 

the corresponding mean-field model. Here the corresponding functions describing the 
conditional expectations can be explicitly computed For more on the analogies be- 
tween non-Gibbsian measures on the lattice and discontinuous behavior of conditional 
probabilities in mean- field models see [Zj and jllj . 

The motivation for using the Morita-approach from the point of view of theoretical 
physics is however to leave these conceptual problems aside and take it as a source 
for approximation schemes [141 115j . Such schemes can be obtained by taking certain 
simplified trial disorder-Hamiltonians that are chosen e.g. demanding that a finite 
number of moments of the distribution of the disorder variables coincide with that of 
the true distribution. Then one would like to solve the resulting Morita approximant 
model and hope that relevant features of the solution are the same as that of the true 
model. 

It might seem hopeless to justify such approximations in general for non-trivial 
lattice models. It is therefore valuable to fully understand at least simple toy models that 
can be explicitly treated. This is want we want to do here. We will give here a complete 
discussion of the quick naive "solution" of the mean- field random field Ising model, 
based on a very simple approximant joint measure containing only one parameter |2.3| . 
This so-called solution is fairly old and the computations are trivial, but the justification 
of the resulting equations is subtle. So it is worth to reconsider it from a rigorous point 
of view and straighten out some wrong claims in the literature, answering a question 
of Kiihn. In particular we take issue with the statement of Ktihn ^3], describing the 
work of [23]. He writes: "it reproduces the exact solution at the cost of introducing a 
single 'chemical potential' to fix the average value of the random field, which creates a 
term in the modified Hamiltonian that introduces no non-locality into the model over 
and above that already contained in the definition of the Curie- Weiss limit. This is a 
remarkable result in the light of concerns raised about the appearance of non-Gibbsian 
measures within the equilibrium ensemble approach and the identification of the RFIM 
as providing a realization of a kind of 'worse-case scenario' in the non-Gibbsian world..." 

We will indeed see that even for this simple model the situation is subtle and the 
validity of the solution is fundamentally related to the analogue of "non-Gibbsianness" 
in the mean- field context. To see this, we will start in Section 2 by reviewing the quick 
"solution", following [Ej. This provides us with two equations for two parameters, the 
magnetisation and the Morita chemical potential. These equation have in fact solutions 
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for which the magnetisation-variable takes the known value of the spontaneous mag- 
netization. However, we note that this solution corresponds to a wrong (metastable) 
saddle-point approximation for the approximant measure and therefore the naive deriva- 
tion given above is flawed. Moreover, we will prove that in the low temperature regime 
it is even strictly impossible to choose a chemical potential such that the magnetic fields 
become symmetric. In brief, the theory based on the Morita approximant measure with 
just one chemical potential fails. 

How can we understand then that the two equations derived by a wrong line of 
argument yield the correct value of the magnetization? Is this just accidental? We will 
see in Section 3 that the same two equations come up in a different way as consistency 
equations for the conditional probabilities of the true joint measures of the model with- 
out approximations. Here however the fixed Morita-chemical potential is replaced by a 
random variable. It is in this context that we will finally understand that the validity 
of these equations and the almost sure discontinuity of the conditional expectations are 
consequences of each other. 



2 Invalidity of single-site Morita approximation approach 
for the Curie Weiss Random Field Ising Model 

We consider the mean-field random field Ising model. It is defined in terms of the 
following formula for the quenched Gibbs expectation for fixed choice of the random 
fields. 

Quenched measure: 

exp(A , (E£=i ff 2 + ^5Z2=i( e, W + h o) <7 ' 

^/3, £ ,h ,iv[ 7 7[l,iV]]( c7 [l,Af]) := 7 f i 

Z /3,e,h ,N[V[l,N]\ 

Here the spins ("dynamical variables") take values <7j = ±1 and the random fields take 
values rji = ±1 with equal probability. We denote their distribution by P. We stress 
that the partition function appearing in the denominator depends on the realization of 
the random fields f][i,N] describing the disorder. We allow from the beginning also an 
external magnetic field ho, but we are mainly interested in the case ho I 0. 

What one understands by the "solution of the model" is the characterization of 
the behavior of this measure on the cr's for a large set of 77's, having asymptotically 
P-measure one. This has been done in great detail |22[IT|l%j. an d so in this model there 
is no need for any approximation based on the Morita-approach in order to solve the 
model. Most basically, we know the phase structure in zero external field, for any choice 
of the parameters (3, e. We recall that for large (3 and small e in zero external magnetic 
field ho the model exhibits a spontaneous magnetization whose value mis a solution of 
the equation 

m = — ^ tanh (3{m + e) + tanh f3{m — e) j (2) 

We also know finer properties of the quenched distribution above, like its dependence 
on the volume label N , for fixed realisation of the random fields. This can be asymp- 
totically described by the in the "metastates formalism" [S], a notion due to Newman 
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and Stein |19l I2U| . For general background on this notion in the theory of disordered 
systems see [2*TT 12*]. 

Knowing the correct solution, our point in this note will be however to put the 
Morita approximation scheme outlined above to the test. Now, in the Morita-approach 
one looks at the joint measures on the product space of the spin variables a and the 
disorder variables rj. They are simply composed from the quenched measures and the a 
priori uniform distribution of the random fields by the following obvious formula. 

True joint measure: 

^/3,e,Ao,iv(^[l,Ar],f7[l,jV]) = V>/3,e,h ,N [V[1,N]} (°"[l,iV]) ( 3 ) 

The approximant measure we want to consider is obtained by putting a single-site 
disorder potential A^ i=1 ?7i with just one free parameter A that has the meaning of a 
chemical potential governing the mean value of the random fields. 



Morita-approximant measure: 

exp 



e x P (m (Eili ^) 2 + P Eili Ni + ha)(Ji + A ^=1 Vi J 

^A;/3,e,ft ,iv(o'[l,Ar]) ? ?[l,Ar]) = 7} (4) 

We stress that the partition function does not depend on Tju^n, in contrast to (^Q). The 
Hamiltonian of this measure contains no non-local couplings of the random fields. 

Then the idea of the naive Morita approximation- approach is as follows: 1) For 
any fixed A, compute the large-iV limit of distribution of this model. 2) Choose A = 
\(P,e,ho) such that the expectation of the random fields coincides with the true joint 
measures, i.e. it vanishes, lim^ J K\.p j£ j l0 ^(dr]i)r]i = 0. More precisely the value of 
A will depend on N, but it will have a well-defined limit as N f oo. 3) Then, the 
distribution of the Morita approximant measure taken with this value of the bias of the 
random fields A, should be close to the true joint measure. E.g. we should have that 
the distribution of the spin average jj Y2i=i a i nas the same infinite volume limit in the 
true joint measure and in the Morita approximant measure. 

Let us write down the following precise formulation in order to have a well-defined 
starting point of discussion. 

Single-Site approximation conjecture: Let (3, e, ho be fixed. Then the conjecture 
is that there is a value A(/3,e,/io) such that lini/v ^A;/3,e ) /i ,A r ( c ^7i) ? 7i = an d that for 
this value we have that 

^ N „ N 



^E^O (5) 



i=i i=i 

Is this conjecture true? How does this relate to the proved a.s. discontinuity of the 
conditional expectations of the true joint measures? Let us review the quick derivation 
of the solution of the model based on this conjecture (we follow here |14|). 
Naive (problematic!) derivation: Look at the partition function of the Morita 
approximant measure, putting ho = from the beginning, and use a simple Gaussian 
identity (Hubbard-Stratonovitch transformation) to write 

*W =2"/ -^^e^(-pm^(rn)) (6) 
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The function appearing in the exponent is TV-independent and is given below in (|24[) 
by putting ho = 0. Using the Laplace method to compute the integral we must have 
■jj^<&\ t /3, e ( m ) = 0. This is an equation for the minimizer m of the form 

E fc=± iSmh(/?(m + gfc))e Afc 
Sfe=±i cosh(/3(m + ek))e Xk 

The parameter A is fixed such that the mean of the magnetic field sum divided by 
N vanishes, i.e. liniTVfoo -§\ log ^A;/3,e,W = 0- This requires at the minimizer m that 
~§\&X;l3,e( r n) = 0. This requires that 

e -2A = cosh(/?(m + e)) 
cosh(/3(m — e)) 

The equation (|SJ) shows that m and A are in one-to-one correspondence to each other. 
From (J2J) and (jSJ) follows the well known (and correct) mean field equation (j2J. So it 
seems that the Morita approximation approach becomes exact in this case and we are 
done. 

Kiihn writes appropriately: This result liffl - simple and reassuring as it is - must 
be regarded as remarkable in the light of concerns raised about the appearance of non- 
Gibbsian measures within the equilibrium ensemble approach and the identi- 
fication of the RFIM as providing a realization of a system exhibiting almost surely 
non-Gibbsian joint measures ffi f77)/. 

Indeed, we note that the "derivation" is flawed because of the following fact. 
Worrisome fact why this derivation is wrong: Suppose that (3 > 1, e > and 
A > are fixed. Then the minimum of the function m i— ► $A;/9,e( m ) is attained at a 
unique positive value m*(A) (as we will see below). Therefore there cannot be a pair 
(m*(A),A) satisfying (|SJ). So the solution (m, A) obtained by (J7J), (jHJ) corresponds to a 
wrong value for the free energy. 

The remaining Morita mystery: Why does the wrong minimizer give the correct 
equation for the magnetization? 

It is the purpose of this note to clarify the situation. We will be even more general 
and more careful here and allow for a possibly non-zero external magnetic field h$. This 
we do in order to investigate whether taking the limit /io | only in the end will help 
us to solve the problem of this approach. 

We can readily solve the model for any choice of the parameters of inverse temper- 
ature /?, strength of random fields e, external field h$ and Morita chemical potential A. 
As usual in mean field models there is convergence to (linear combinations) of product 
measures over the sites i. Indeed, any limit measure must be a mixture of product mea- 
sures. This is clear by de Finetti's theorem since the limit of exchangeable measures 
inherits the property of exchangeability. 

Now, solving our simple model is almost trivial when we note that by summing over 
the rj first we obtain a resulting effective Curie- Weiss Ising model with a new effective 
homogenous magnetic field acting on the cr's. The computations are simple and will be 
given below for the sake of completeness. Before we do so let us however state the most 
important consequence of this in the present context. 

Theorem 2.1 (Impossibility of single-site approximation of true joint mea- 
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sure) Assume that (3 > 1 and e > are fixed. Then 

{h Gi,3AGK: lip Kx;f3, £ , ho , N ( d v) = 

Af|oo ) ^ 

-a(/3,e),+a(P,£ 

where a(/3, e) is strictly bigger than zero. Here the symbol lim denotes a weak limit. 

In words the theorem states that the set of external homogenous magnetic fields ho for 
which there exists a "compensating" Morita-field A that reproduces the neutral i.i.d. 
distribution for the random fields is bounded away from zero for any (3 > 1. This means 
that the approximation scheme must necessarily fail in the relevant low temperature 
regime: First of all, in zero external field ho it is impossible to produce asymptotically 
symmetric i.i.d. random fields by an appropriate choice of A. This result however, might 
not be too surprising. But the theorem says more: Even choosing ho strictly positive 
and letting it tend to zero afterwards won't help us. 

Having said this, it is interesting to investigate the set of parameters for which the 
distribution of random fields becomes neutral i.i.d. in more detail. Let us make the 
following definition. 

Neutral Set: Fix the inverse temperature j3 > and e > 0. We call the parameter 
set 

K(j3,e) := {(h ,X) «x R, Urn K x . AeM;N (d V ) = P(cfy)} (10) 

the neutral set. Obviously TZ(f3,e) = —1Z{(3,e) by the symmetry of the model. 
Then we have the following theorem. 

Theorem 2.2 (Structure of neutral set) Assume that (3 > and e > are fixed. 
Then the set TZ((3, e) is the union of two semi-infinite curves, related to each other by 
reflection at the origin. They are connected if and only if (3 < 1 . 

More precisely these curves are of the following form. 

There is a continuous increasing function ho i— > lp^ £ {ho) that is defined an open 
interval of the form (a((3, e),oo) and takes positive values. The left endpoint of the 

interval satisfies a((3, e) \ > ^ ^ > 

1 = for (3 < 1 

Define 

TZ + (P,e) := {(h ,-l^ £ (ho))\ho & (a(/3,e),oo)} (11) 

Then, the neutral set has the form 

(K+(J3,e)\j(-K+(J3,ej) for (3 > 1 

[K+((3,e) U (-ft+GM) U {(0,0)} for f3 < l 

Hence, for [3 > 1 the neutral set is disconnected. For (3 < 1 we have moreover 
l*p e (0+) = 0, and hence the neutral set is connected. 
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This result on the neutral set is a consequence of the solution of the Morita approxi- 
mant for any choice of the parameters. We will now describe the behavior of the Morita 
approximant for general choice of the parameters. Then we will derive as a conclusion 
the explicit condition for the neutral set. 

We need some definitions. Define the effective magnetic field-like parameter 

h = h Q + hp, e {\) (13) 

with the function 

7 tx\ 1 i cosh(A + /3e) 

h,sW := ~2f3 g cosh(A - (3e) (U) 

Define the joint single-site measures depending on the parameter set, and on an 
additional (magnetization-like) parameter m € M. 

exp(/?((m + er]i + h )(Ji + Xrji)^ 

^, h M^m) - 2Efc=±iC08h( ^ (m + efc + ^ ))c ^ (i5) 

_ exp(/3(m + h)ai) exp((/3eo-j + X)rji) 
2 cosh /3(m + h) 2 cosh(/3ea i + A) 

Here we found it convenient to express the joint distribution on the r.h.s. appearing 
under the i-product in the form Prob(cij, rji) = Prob(o"j)Prob(r7j|cjj). In this way the 
marginal on the cr's can be readily read off. We see that the role of the parameter h is 
to provide an "effective magnetic field" acting on the spins. 

Theorem 2.3 (Solution of Morita approximant) Assume that the parameters 
(3, e € (0, oo) and A, ho € (— oo, oo) are fixed. 

(i) : Assume at first that h ^ 0. Then we have the weak convergence 

N 

1™ K\;f3,e,ho,N(v[i,No],V[i,N ]) = Yl ^Weto [ mCW (Pi h )] Vi) ( 16 ) 
|oc i=i 

Here we have denoted by m cw {j3,h) the solution of m = tanh(/3(m + h)) that has the 
sign of h, for h ^ 0. 

(ii) : For h = we have the weak convergence to the symmetric linear combination of 
product measures 

^^X;p,e,hoM a [l,NohV[l,N }) 

, No N (17) 

= 2 II ^;P,eM [m CW {^ 0+)] (a t , Vt ) + - JJ n x . Mo [m cw (P, 0-)] (a t , m ) 
i=i i=i 

Remark: Of course m cw (/3, h) is the magnetization of an ordinary Curie Weiss Ising 
model in an external field. 
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Proof: Let us write the Morita approximant joint measure as a marginal on the it's 
times the conditional measure of the random fields given the cr's, that is 

K\;f3,e,h ,N( a [1,7V] ) V[1,N] ) 

ex p( 22V (Sill a if + P h o YhLi a % + Ya=\ lo g cosh(/3eo-i + A) J JL e (Pea+X)m 

Norm . AJ- 2 cosh^eo; + A) 

(18) 

This shows us that the marginal distribution on the cr's is given by an ordinary ordered 
mean field Ising model of the form oc exp (SfcLi a i) 2 + (3h YliLi a *) /Norm . with the 

effective field h. From here the limit statements are obvious by the known convergence 
results of the Curie Weiss Ising model to the corresponding (linear combination of) 
product measures. □ 

Now, from the explicit solution we may derive explicit information on the neutral 
set. In order to do so, note at first the elementary properties 



| — e, for A I -co 

V.( A ) = { 0. forA = 
| e forA | oo , 

and it is a monotonically increasing in A and odd. It maps R to the interval (— £, e). 

Theorem 2.4 (Explicit description of neutral set) Assume that (3 > and e > 
are fixed. The decomposition hlty) holds with 



m cw ([3M-hAl)) 



sinh(2Z) i ( 19 ) 



sinh(2/3e) 



This set can be written as a graph in the form with a continuous increasing lp e {ho) 
that maps the interval (a(/3, e),oo) onto the interval (hp e (a(f3,e)),/3e) where a(/3, e) is 

uniquely given by hp je (a(/3,ej) = \ sinh -1 ( sinh(2/?e) m cw (f3, 0+) J . 



Remark: Note that the above expression for a(j3, e) implies that a(/3, e) = if and 
only if the spontaneous magnetization m cw (/3, 0+) vanishes, i.e. /3 < 1. 

Proof: Suppose that j3 > 1. Then, in order to have convergence to a symmetric product 
measure on the random fields we must have that the parameters A; 0, e, h^ are such that 
fe / 0. (Indeed, for h = the distribution on the spins converges to a symmetric 
mixture of two different product measures. But from this it is obvious that also the 
random field distribution will be a mixture between two different product measures.) 
This shows that (0, 0) 7Z(P, e) in that case. 

Suppose however (3 < 1. Then (ho, A) = (0,0) implies h = which implies that the 
distribution of the cr's is a symmetric product measure. But this implies that the distri- 
bution of the random fields will be a symmetric product measure so that (0, 0) 7£(/3, e) 
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in that case. 



So, we are left with the case h ^ 0. We can treat the cases (3 > 1 and (3 < 1 on a 
unified basis. Now, conditional on the value of a the 77, have an expectation value of 
tanh(A + (3eo~i). We use the simple identity 

f h^-ufl ^ -^(1 ~ L 2 )ai + L(l — B 2 ) 
tanh(A + /fed;) = 1 - ff2L 2 where 

L = tanhA, 1? = tanh/3e 

for cjj = ±1. So the distribution on the random fields rji converges weakly to a product 
measure with individual expectation value 

V Ik <a \ B(l-L 2 )m™((3,h) + L(l-B 2 ) 

km J Kx;f},e,h M d Vl)Vl = 1 _ B 2 L 2 ^ 

Put I = — A and use tanh(/)/(l - tanh 2 (/)) = sinh(2/). So, in order to have the 
desired convergence to the symmetric product measure we must have 

This equation can only hold if ho — hp >£ (l) and I have the same sign. By symmetry we 
can assume that I > 0. But this implies that ho > ( since hp jS (l) > 0.) 

So, it suffices to look for all pairs (I, ho) with I > that satisfy the consistency 
equations (|22|) . The small trick we are using now is to fix the I and ask for ho rather than 
doing it the opposite way. Fixing I we see that the l.h.s. runs monotonically through 
the open interval (m cw (/3, 0+), 1) when ho runs in the "allowed range" (ha e (Z), oo). 

So, the set of I > such that there exists a solution ho is determined by the condition 
(m cw (/3, 0+), 1) 3 S in^2^e) • Equivalently, this is the open interval I € (a(/3, e), /3e). 
Moreover the map to ho is continuous and monotone by known properties of the function 
m (P, h). So it can be inverted and this yields the claim. □ 



So what has happened in the naive (but wrong) derivation of the mean-field equa- 
tions (JJJ) and l(8|)? In order to see this let us write down a representation of the finite- N 
approximant measures. As a result of a Gaussian transformation on the level of mea- 
sures we get the following formula. 

Proposition 2.5 In finite volume N we have the identity 

I dmex.p(—/3N& j 3 e f l() (in)) 

K\;i3,e,h ,N(<r[l,N],V[l,N)) = T7~~ 1 a^r^ TTTT I I ^X;/3,e,ho l m \ v^ii Vi) ( 23 ) 

J dmexp(-(3N$x;f3,e,h { m )) t=i 



Here 

2 [3 



®X;f3,e,ho( m ) = ~7T ~ \ lo S Yl cosh (^( m + ek + h o))e Xk 



m 2 1 



2 p 

where Const (/3,s) does not depend on m. 



k=±l 

log cosh (/3(m + h)) + Const (/?, e) 



(24) 
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Remark: The second equality for <3> can be seen e.g. by reexpressing the first cosh as 
a sum over a spin s = ±1 and exchanging the s and A;-sums. 

Proof: We use a Gaussian transition kernel from the cr-variables to an auxiliary real 
valued variable m given by T(dm|(T[ lj jv]) = exp^— ^j-(rn — a% \ ^<im/Norm.. We 
define a "big joint measure" on the spins, the random fields and also the auxiliary 
magnetization-like continuous variable by the formula 

Mx^,e,h ,N{dm,a[ 1>N ],r][ 1>N ]) := Kx-/3,e,h ,N{^[i,N],V[i,N])T(dm\a[i tN ]) (25) 

We see that m concentrates very nicely around the value of the empirical average of the 
true spins in this measure. Then the non-normalized density of this "big joint measure" 
is given by exp^— @j-m 2 + f3 J2i {( m + er li + ho)ai + Xtj^. Use this to express the "big 
joint measure" in the form of a marginal on the m times a conditional measure on the 
(a, rj) given the m. Prom here it is simple to get the desired formula. □ 

So, conditional on a value of m, the pairs (<7j, rji) are independent. We have then for 
their conditional mean values 



Z k=±lS mh(f3(m + ek + h ))e Xk 
a W0H(ffl)ffl = E^cosh^ + ^ + ^e- 
v—v r , Efc=±i^cosh(/3(m + eA: + / i o))e Afc 



(26) 



We remark that, with this notation, we have that (the version for general ho of) the 
saddle point equation Q is equivalent to the consistency equation for the magnetization 
written as 



m 



CTl=± 

The (version for general ho of) the neutrality equation Q is written as is equivalent to 

(28) 



= Y 7r A;/3,£,/ l0 H('7l) ? 7l 
771 =± 



Now, the large- limit of the model is obtained by looking at the absolute minimizer 
of the function m 1— > <3?(m). But note that the representation for $(m) given in the 
second line shows that is has the double-well form of the corresponding function in an 
Ising model in the external field h. It is an elementary property of this function that 
its absolute minimizer has the same sign as h. But this shows that the relation (jHJ) 
can not be true for the absolute minimizer. Instead the solution of 0, (jHJ) corresponds 
to the second local minimum which is not the absolute minimum but the metastable 
minimum. 
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3 Validity of consistency equations and almost sure dis- 
continuity of conditional expectations 

So how can we understand the fact that the correct solution of the model is obtained 
by solving equations (|27|) and (|28(1 although the solution corresponds to the wrong 
saddle point? The solution to this puzzle is due to the fact that the naive equations 
have rigorous counterparts in the following sense. The equations we are going to state 
now appear as consistency equations for the conditional probabilities of the true joint 
measures. 

Proposition 3.1 (Consistency equations for true joint measure) There is a 
function ^N(V[2,N])t depending on the parameters 0,e,ho, which is invariant under per- 
mutation of (rji)i = 2 t ... n such that we have 



y~i-K/3,g,feo,iy(°"i) (7 i 

CT1 

K P,e,h ,N{o- [2 ,N],Vl2,N]) ( ^X N (rj [2iN] );/3,e,ho [jj 

T [2,N],V[2,N] \CTi=± 1=2 / 



(29) 



N 



0= S K ^,h ,N(o- [2 ,N],ril2,N]) ( Yl ^(^JsiS.e./io^Z)^]^)^) (30) 
a %N] M[2,N] \?7i=±l i=2 J 

Proof of the proposition: The proof is based on the following lemma. 

Lemma 3.2 (Representation of conditional probability of true joint measure) 

The single-site conditional probabilities can be written in the form 

1 N 

K{3,e,h ,N (0-1 , 771 |<7[ 2 ,AT] , ?7[2,JV] ) = 7T Ajv ( V%N] );/3,e,h [jj Y a i\ 1 Vl ) ( 31 ) 

i=2 

where 

w x 1, Z P,e,h MVl = ~,V[2,N]} . . 

Ajv IV[2.N] = o !og -z f : ] 32 

2 Z^ eM ,N\r)\ = +,V[2,N}\ 

Proof of the lemma: By a simple computation we have for the single-site distribution 

K/3,e,h ,N {o-\ , Tji \o [ 2 ,N] , V[2,N] ) 

1 ( a( 1 A , ,,x 1. Z /3,e,hoMVl = ~,V[2,N]] \ ( 33 ) 

= Tf exp /?(— > +£7/i +/l J<7i + -log- f ■ r X Tjl J 

Norm V v iV ^ 2 ZfrsMMVl = +, ?7[2,JV]] y 

and this shows the claim. □ 
Continuing with the proof of the proposition we use the formula for the conditional 
probabilities writing 



1 N 

Y K P,e,hoM a [2,N},V[2,N]) ^\ N (v%N]);^M [jy= X) aj \ ( CTl > 

CT [2,JV]i'7[2,]V] «=2 



(34) 
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But this equation gives the equation for the magnetization (|29f) by summing over a±. 
Using the symmetry of the distribution of r/i we get (|3T)|) . □ 



Let us now summarize what we know by the rigorous solution of the random field 
model about the limiting distribution of the pair of random quantities entering the 
single-site kernel tt. In words, the distribution becomes sharp in the case of non-zero 
external field. It becomes sharp but double valued in the case of vanishing external 
field. In view of the last lemma this statement is a different way of saying that there 
is a jump of the conditional probabilities when the empirical random field sum of the 
conditioning is infinitesimally perturbed around its typical value 0. Now, the rigorous 
statement is as follows. 



Theorem 3.3 (Convergence of true joint measures) 

(i) Suppose that ho > 0. Then we have the weak limit 



Here (m* (ho) , X* (ho)) is a solution of the consistency equations \21\j and i28\) . 
(ii) Suppose that ho = 0. Then 



hmi^ 0=0 ,. 



1 N N \ 

Y J2 ai G " ' X n(V[2,n]) € • ^Vi > 

i=2 i=l J 



x <5 A 



where (m*,X*) is the unique solution of the consistency equations \2% and H28\) with 
m* > (and, as a consequence X* < 0). 
As a consequence we have 

fl N \ 1 1 

NpL K(3 ' e ' ho=0 ' N ( N G ' ' Xn ^1 2 ' N ^ £ J 2 6m * X 6x * + 2 5 ~ m * X 5 ~ X * 

Remark: We see that the system chooses the particular value of Aat^^tv]) (that has 
the opposite sign of the magnetisation) itself! 

Proof: We only sketch the proof. We rewrite the quotient of partition functions 
appearing in the definition of Ajv^^at]) in the form 

X N{V12,N}) = 2 log/^^jvfoi = +,7/[2,7V]](exp(-2/3ecri)) (35) 

From here Theorem 3.3 follows from the work done for the quenched model in [Sl^J. Let 
us focus here only on the interesting case of vanishing external magnetic field ho = 0. In 
this case it was shown that, under the condition of positive sum of the random fields the 
empirical average of the spins concentrates sharply around the positive magnetisation 
m* (positive solution of ((2j)) w.r.t. to the quenched Gibbs probability. (This is true 
for "typical values" of the random field sum, that is for N^ s < Vi < N ^ + \ and 
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these values get all mass w.r.t. P in the large- N limit). At the same time the quenched 
Gibbs probability /X/3, e ,/i ,Ar[??i = +, V[2,N]\( cr i = +) aquires a sharp value that is related 
in a simple way to m*. From Q35|) this gives the value of A*. □ 

Not assuming the knowledge of the solution of the quenched model we can reverse the 
argument in the following way in order to solve the model. Look at the consistency equa- 
tions for the true joint measure (|29 |) .(|3l)| ) . Take ho > 0. Then it is very plausible with- 
out much a priori knowledge that the distribution of the pair Yli!=2 °»> ^n(i][2,N])) 
under the true joint measure should converge to a Dirac measure S mj \. (This is in par- 
ticular clear, if we assume that A at has the form Q35|) and assume that the quenched 
magnetization becomes sharp for typical realization of the random fields in a positive 
homogeneous external field.) But this means that the outer integrals in the rigorous 
consistency equations become sharp. So, the limiting value (m, A) must necessarily sat- 
isfy the naive consistency equations (|27|) and (|28|) . These equations can then be solved 
and afterwards we let the external magnetic field /io tend to zero from above to discover 
the known solutions for the model. 

Let us finally see that, in the case of ho = the validity of the naive equations implies 
that there must be discontinuous behavior of the conditional expectations as a function 
of the average of the random fields appearing in the conditioning. Indeed, suppose that 
Xn(i][2,n]) were a continuous function of 2^^2 r ?«- Then, by the law of large numbers 
it would have to be constant in the large- N limit. But by reasons of symmetry this 
constant would have to be zero in the case of hg = 0. But this is in contradiction to 
the non-trivial solution of the naive equations and (|5j). To summarize the last line 
of argument in catchy terms: Non-Gibbsianness is necessary to help the metastable 
solution to provide the right answer. 

We remark that the purpose of this note is not to attack the Morita approach in 
general as a valuable heuristic method in theoretical physics to predict the behavior of 
disordered systems when a rigorous analysis is not available or not yet available. 

As pointed out to us by Reimer Kiihn, one could also argue that the second of the 
naive equations (jSJ), which demands that m and A at the physical fixed point must 
have opposite sign, renders the region of integration for the partition function © which 
includes the other fixed point as unphysical and so to be excluded from the domain 
of integration. This line of reasoning would render the "naive" argument correct and 
this is not the first occasion in physics where such things happen. While there seems 
no direct rigorous justification for this procedure we have shown that one is able to 
understand the validity of the naive equations by viewing the parameter A properly 
as a stochastic quantity. This might give hope that results obtained by approximation 
schemes based on the Morita approach provide correct answers also in more complicated 
situations where a rigorous analysis is lacking. A better understanding of this would 
pose a fascinating challenge. 
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